Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



A Measurement of the Ratio of the W + 1 Jet to Z + 1 Jet Cross Sections with 

ATLAS 

A. Meade and B. Brau on behalf of the ATLAS Collaboration 
Department of Physics, University of Massachusetts, Amherst, MA, USA 

The measurement of hadronic activity recoiling against W and Z vector bosons provides an 
important test of perturbative QCD, as well as a method of searching for new physics in a model 
independent fashion. We present a study of the cross-section ratio for the production of W and Z 
gauge bosons in association with exactly one jet _Rjot = a(W + l]et) /a{Z + Ijet), in pp collisions at 
y/s = 7 TeV. The study is performed in the electron and muon channels with data collected with the 
ATLAS detector at the LHC. The ratio Rjct is studied as a function of the cumulative transverse 
momentum distribution of the jet. This result can be compared to NLO pQCD calculations and 
the prediction from LO matrix element + parton shower generators. 

I. INTRODUCTION 

Measurements of the cross sections of the W and Z bosons in association with hadronic activity are important 
tests of the standard model, particularly perturbative quantum chromodynamics. In addition many searches for 
models of physics beyond the Standard Model have significant LF+jets or Z+jets backgrounds, adding to the 
importance of understanding these processes. Individual measurements oiW oi Z cross sections in association 
with jets are limited by systematics that are shared between the two measurements, such as the luminosity and 
the jet energy scale. In contrast, in measurements of the ratio between W and Z cross sections, many of these 
uncertainties cancel due to the similar nature of W and Z production. For this reason the W/Z cross section 
ratio can be the basis for a program of high precision measurements. 

This study reports on a measurement of the ratio of production cross sections of the W and Z bosons with 
exactly one associated jet with > 30 GeV as a function of jet transverse momentum threshold. This binning 
implies that each bin contains all events with a jet above the threshold. The data for this analysis consists of 
the entire 2010 ATLAS dataset, operating at the Large Hadron Collider (LHC). This corresponds to 33 pb~^ 
of integrated luminosity, taken at a center of mass energy of 7 TeV. A paper on this measurement has been 
produced by the ATLAS collaboration [P. 

The ratio is measured in a fiducial measurement volume of the detector where the leptons and jets are 
measured with good resolution, in order to avoid theoretical uncertainties associated with modeling particles 
outside of the detector volume. The fiducial region is defined by the following kinematic ranges of the lepton(s) 
I, neutrino v, and jet: p^r^ > 20 GeV, > 25 GeV, p^^' > 30 GcV, Irf"^] < 2.8. The electrons and muons 
were restricted to different pseudorapidity ranges, for electrons 1.52 < \ri\ < 2.47 or \ri\ < 1.37, while for muons 

I77I < 2.4. For the W, the transverse mass, defined as mx = y^2p;^p^(l — cos((/)^ — (j)'^)) was required to satisfy 
mx > 40 GeV. The dilepton invariant mass of the Z was required to be within the range 71 < m« < 111 GeV. 

Particle level jets were defined as jets reconstructed in simulated events by applying the anti-fct jet recon- 
struction algorithm [2] with a radius parameter R — 0.4 to all final state particles with a lifetime longer than 
10 ps (including muons and non- interacting particles). Particle level electrons were defined by including the 
energy of all radiated photons within a cone of AR = 0.1 around each electron (with AR = ^ (Ayy)^ -I- (A</))^). 
In the muon case the "bare" final state muon is taken as the particle level muon, with no radiated photons 
included. 



II. THE ATLAS DETECTOR 

The ATLAS detector is described in detail elsewhere [5- Immediately surrounding the interaction point is 
the inner detector (ID), providing precision tracking and vertexing capabilities, surrounded by a large solenoidal 
magnet. Outside of this volume lies the electromagnetic and hadronic calorimeters, used for energy measurement 
and electron and photon identification. The outermost portion of the detector is the muon spectrometer (MS), 
based on three large superconducting toroids and a system of three stations of trigger chambers and precision 
tracking chambers. 
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III. SIMULATED EVENT SAMPLES 

Simulated event samples were used to correct signal yields for detector effects, for some of the background 
estimates, and for comparison of the results to theoretical expectations. Samples oi W ^ iv + A'parton and 
Z ^ ££ + A^paiton (where £ — e,/i, r) were generated using Alpgen [4 . Background and additional signal 
samples were generated with Pythia 5^ and PowHEG [3]. To take into account overlapping pp collisions some 
samples were generated with multiple non-difFractive scattering events overlapping with the primary collision 
event. These samples were reweighted such that the distribution of the number of reconstructed primary vertices 
matched that of the data sample. GEANT4 [7] was used to simulate the detector response for all samples, and 
these samples were subject to the same reconstruction and analysis chain as the pp collision data. 

Predictions for the W and Z cross sections as a function of jet pt threshold at NLO were obtained using 
MCFM, with corrections to particle level calculated using Pythia, to account for initial and final state radiation, 
hadronization, and underlying event. 

IV. DATA AND EVENT SELECTION 

Kinematic requirements for reconstructed muon and electron candidates match those of the fiducial definition: 
Pt > 25 GeV, \ri\ < 2.4. In the electron channel events were triggered based on the presence of an electromag- 
netic cluster with Et > 15 GeV. Electron candidates were required to satisfy lateral shower containment, shape 
and width criteria, and minimal leakage into the hadronic calorimeter in order to be classified as "medium" . 
Consistency between track pt- and cluster energy as well as tighter hit requirements are additionally required 
to categorize an electron as "tight". For the electron channel the pseudorapidity range region 1.37 < jryj < 1.51 
is rejected due to a gap in the calorimeter coverage. 

In the muon channel, events were selected by the trigger based on muon spectrometer hits consistent with a 
track of 10 or 13 GeV, the tighter requirement applying to later data with higher instantaneous luminosity. Muon 
candidates were required to have inner detector and muon spectrometer tracks and pass kinematic requirements 
to be classified as "medium". To be classified as "tight," consistency is demanded between ID and MS track, 
and requirements on the number of inner detector hits were made. In addition, the impact parameters of the 
muon are required to be consistent with the interaction point, and the inner detector track is required to be 
isolated from other energetic tracks, with J^Pt of all tracks within AR < 0.2 less than 1.8 GeV. 

Each event was required to have a reconstructed primary vertex with three or more tracks, and to be within 
150 mm of the center of the detector along the beamline. W candidates were required to have exactly one 
"tight" lepton with no additional "medium" leptons, E^^^^ > 25 GeV, and AIt > 40 GeV. Z candidates are 
required to have exactly two leptons of opposite charge, one "tight" , and one of "medium" or higher quality. 
The invariant mass mu of the lepton pair is required to be in the range 71 < mu < 111 GeV. 

Reconstructed jets were defined using the anti-Zct jet reconstruction algorithm with a radius parameter R — 0.4 
as in the fiducial definition. These jets were required to be within the kinematic range px > 30 GeV, |r/| < 2.8. 
Events were rejected if identified as likely to contain a jet from detector noise or a cosmic ray. Events with jets 
within a cone of AR < 0.6 from an electron are rejected to avoid distortion of the shower shape from the nearby 
electromagnetic shower, which could result in reconstruction inefficiency or incorrect energy measurement. Jets 
from additional interactions in a bunch crossing were rejected by requiring that 75% of the scalar sum px of all 
tracks associated with a jet came from tracks originating from the same primary vertex. This approach is termed 
the jet vertex fraction (JVF) algorithm. Only events containing exactly one jet passing these requirements were 
selected as candidate events. 

Using these W{Z) selection criteria, 12112 (948) events and 12995 (1376) events were found in the electron 
and muon channels respectively. 

V. BACKGROUND ESTIMATION 

Backgrounds to the measurement are determined independently for each bin of the measurement in jet pt 
threshold. These backgrounds are categorized as resulting from electroweak or QCD multijet processes, and are 
subtracted from the total number of events passing selection using the following formula to obtain the signal 
event yield N^ig: 



— Ntot • (1 — /multijct)(l — /cwk) 



(1) 
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The predicted number of events passing selection cuts is shown for all signal and background processes in 
Table H 

The electroweak background fraction /ewk is estimated from simulated data as a fraction of the total multijet 
subtracted event count. This has the advantage of having no dependence on the luminosity, and reduces 
dependence of the background on our measurement volume, as the electroweak processes have similar kinematic 
distributions. 

The multijet background processes are estimated using data driven methods that vary based on the channel. In 
the W ^ ev channel, signal and background templates were fit to the control region 15 GeV < E^^^^ < 55 GeV 
to determine the background fraction. The signal and electroweak background template shapes are taken 
from simulation. The multijet event shapes are found from data by selecting "medium" electron candidates 
and requiring that 2 of the "tight" selection requirements, as well as isolation, are failed, producing what is 
essentially an electron-like multijet sample. The Z ^ ee background estimate was performed using a similar 
method by fitting templates to the dilepton invariant mass spectrum. Because of smaller background, a looser 
electron selection and reversal of two of the "medium" selection criteria were used in the determination of the 
multijet template. 

In the W ~> fiv channel, the multijet background is estimated by measuring the number of events passing 
all selection, and all selection but isolation. The efficiency of this selection is then measured on data control 
samples in the data. Dimuon events from Z events are used to measure the efficiency for muons from electroweak 
processes like the W and Z signal, and dijet events are used to measure the efficiency for muons from the 
multijet background. Since these two efficiencies are significantly different, the multijet background fraction 
can be determined. In the Z iifi channel, the multijet background is very small. In this case a scale factor 
is derived by comparing non-isolated dimuon pairs in data and simulation, and used to normalize a simulated 
multijet sample. 

In the electron W channel, multijet backgrounds arc the dominant background due to the relatively high 
probability of a jet fragmentation product being reconstructed as an electron. The electron multijet background 
fraction was about 16% for the lowest jet threshold, 30 GeV. Z ee, W ti/, and tt events are all significant 
sources of electroweak background, which is around 3.4%. The dominant backgrounds in the W ^ channel 
are W ^ tv decaying to a muon (2%), Z ^ with one muon failing to be reconstructed (3%), and multijet 
events (3%). In both the electron and muon Z channels, backgrounds were very small, less than 1% for all Z 
candidates. 



Process 


W ^ ev 


Z ee 


Process 


W ^ liv 


Z — >■ ^/x 


W eu 


9340 ± 40 


3± 1 


W ^iv 


11860 ±40 


4± 2 


Z ^ ee 


106 ±3 


880 ± 10 


Z — )■ 


360 ±6 


1370 ± 40 


W ^ TV 


191 ± 6 


0.2 ±0.2 


W ^ TV 


234 ± 6 


0.3 ±0.6 


tt 


33 ± 1 


1.9 ±0.2 


tt 


35 ± 1 


3± 2 


Z ^TT 


19 ± 1 


0.3 ±0.1 


Z ^TT 


22 ± 1 


0.3 ±0.6 


Multijet 


1800 ± 60 


2.9 ±0.6 


Multijet 


380 ± 70 


4±4 


Other 


58 ±2 


1.1±0.1 


Other 


117± 1 


8±3 


Total 


11550 ± 70 


880 ± 10 


Total 


13010 ±80 


1380 ± 40 


Data A^tot 


12112 


948 


Data A^tot 


12995 


1376 



TABLE I: Predicted and observed events in data in the electron and muon channels for the W and Z selections for 
J Cdt = 33 pb~^ . Background estimates are quoted for a jet pr threshold of 30 GeV. Only statistically errors are 
displayed. "Other" includes contributions from diboson and single top events. 



VI. CORRECTION PROCEDURE 



Event yields were corrected for trigger efficiency (e^j.;^), lepton identification efficiency (e ), and differences 
between boson acceptance at detector level (Aioco) and particle level (^part) due to resolution effects (Cy = 
). The number of signal events for each boson at particle level was then obtained using 

jy'^-y = ^ (2) 

^'part J V V r'^ ■ 
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FIG. 1: Detector jet spectrum correction (Cj^^) on _Rjct in the muon channel derived from Alpgen. The uncertainty, 
shown as a dashed red hne accounts for the difference between Pythia and Alpgen generators. The black error bars 
show the statistical uncertainty alone. 



where the boson corrections Cy correct the observed phase space to the truth (particle) level kinematic phase 
space, accounting for the resolution of leptons and E^^^^. 

Trigger and identification efficiency for leptons were determined from an unbiased control sample by selecting 
a well identified "tag" lepton and a additional "probe" inner detector track in Z ^ ££ candidate events. 
The efficiency is the fraction of probe tracks matched to a reconstructed lepton and passing all identification 
requirements. These efficiencies were found to be largely independent of jet kinematics, other than through 
correlations with lepton kinematics due to the jet recoiling against the W or Z. Therefore, these efficiencies 
were binned in lepton and rj for the electron channel, and 77 and (f> for the muon channel, for all jet px 
thresholds. 

When the ratio of cross sections is measured, jet resolution and scale effects almost completely cancel. A 
correction C^g^ was applied to take into account small differences. This correction is defined as the ratio between 
the i?jot value as a function of particle level jet threshold and the i?jot value as a function of reconstructed 
jet pt threshold. This factor is measured in simulation samples and applied as follows: 



-^jet = — ^ X C^^^. (3) 

part 

This correction is shown for the muon channel in Figure [T] Differences from unity correspond to an offset in 
the amount of migration between jet pt bins between W+jet and Z+jet events. This is mostly due to selection 
differences in the two channels before the jet selection, such as the requirement that the decay products be 
located in the detector volume. Note that binning in terms of jet threshold eliminates migration across the 
upper bin edges. 



VII. SYSTEMATIC UNCERTAINTIES 



In order to properly take into account correlations between the W and Z measurements, and to maximize 
cancellation of systematic uncertainties, all such uncertainties were measured as a relative change in i?jet- 
The total systematic uncertainty varies from 4% to 15%, generally increasing with increased jet threshold. 
Systematic uncertainties on the measured ratio were divided into five broad groups: multijet background, 
electoweak background, boson reconstruction (combining e^^.^^ x x Cy), the jet spectrum correction C'j^tJ ^^'^ 
uncertainties related to generator differences, as shown in Figure [2] 

Systematics on the multijet background were estimated by varying the parmeters of each method, such as the 
selection of control samples used to derive the background fractions. Systematics on the electoweak background 
are estimated by varying the lepton identification criteria, by varying momentum scale and resolution, and by 
comparing samples generated with and without additional overlapping pp interactions. 
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FIG. 2: Relative systematic uncertainties on Rj^t in the electron channel (left) and in the muon channel (right). The 
top plot displays the total systematic and statistical uncertainty (shown as red dashed line) versus jet px threshold. 
The lower plot shows the breakdown of the systematic uncertainties. The "boson reconstruction" entry contains the 
uncertainties related to the leptons and _E™'°'^ (including trigger and lepton identification). The "jet" entry contains 
systematics of the jet correction as well as the jet energy scale and resolution. Uncertainties from each group were added 
in quadrature. 

Boson reconstruction systematic uncertainties included those associated with identification and trigger effi- 
ciency, as well as scale and resolution effects. The single lepton trigger and identification efficiency uncertainties 
were derived using the methods documented in Ref. [5], and then propagated to uncertainties on i?jet- This 
leads to a total identification efficiency uncertainty of 1.1%(1.7%) for electrons (muons) independent of jet px- 
Trigger efficiency uncertainties were very small in both channels. 

Uncertainties involving the scale of lepton momenta and E^^^'^ were evaluated by directly scaling up and down 
the quantity in question by the measured scale uncertainty. In each case the deviation of i?jet was applied as 
a systematic uncertainty. Similarly, resolution effects were evaluated by smearing the quantities by a gaussian 
with the width of the measured resolution, repeating this process and taking the root mean squared deviation 
as the systematic uncertainty. Lepton momentum scale and resolution measurements and their uncertainties 
are determined by comparing the width and position of the invariant mass spectrum of Z candidates in data to 
that in a simulated sample. 

Jet related systematics included those involved with jet energy scale (JES), jet energy resolution (JER) and 
their associated uncertainties. These quantities were determined from data and simulation comparisons [S]. 
The JES uncertainty includes components from calibration and jet sample composition differences. The JES 
calibration uncertainty varies with |?7| and px, and ranges from 4% to 8%. The JES and JER were measured 
with di-jet events, which have different proportions of quark and gluon initiated jets than events containing 
vector bosons. Therefore, an uncertainty was assigned to account for the difference in calorimeter response 
between jets in y + jet events and the di-jet events used for calibration, ranging from 2 to 5%, and was added 
in quadrature to the JES calibration uncertainty. The total JES uncertainty ranges from approximately 10% at 
20 GeVto 5% at 100 GeV. The JER is measured to vary from 2% at 20 GeVto 1% at 100 GeV. The uncertainties 
on i?jct due to JER and JES were evaluated using the same approach as for the lepton scale and resolution 
uncertainties. The uncertainties on Rj^t due to the JER and JES were found to be approximately 0.5% and 2% 
respective. This uncertainty includes a very small term to to take into account the uncertainty on JER itself. 

Jet related systematics also included the dominant uncertainty related to multiple pp collisions, the uncer- 
tainty of the efficiency of the JVF algorithm. Applying this algorithm to simulation with multiple interactions 
was found to produce consistent results with simulation not including additional interactions. The residual 
difference on i?jet between these cases was used as a systematic uncertainty on C^^^. 

To account for systematics associated with generator modeling, correction factors were computed with samples 
generated with Pythia instead of Alpgen, and the observed difference was applied as a systematic uncertainty. 
Systematic uncertainties were assigned from this variation to the following corrections: (Cj^et)j ^^e boson re- 
construction correction Cy, and the electroweak background estimation /cwk- At large jet px threshold, where 
the statistical uncertainty on the measurement dominates the total uncertainty, this systematic uncertainty is 
limited in statistical precision due to the size of the samples used, and is the dominant systematic uncertainty. 
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Systematic uncertainties were also estimated on the theoretical prediction. As our NLO parton-level cal- 
culation does not include the effects of hadronisation and underlying event, a correction was computed using 
Pythia as a function of jet threshold. The uncertainty on this correction was evaluated by observing the 
effect of various generator tunes on the final i?jot result. These tunes increased or decreased the amount of 
underlying event, or varied the parameters controlling initial and final state radiation. Renormalization and 
factorization scale uncertainties were also included in our systematic uncertainties, as were PDF uncertainties. 



VIII. RESULTS 



The ratio i?jct was measured in the fiducial region of the ATLAS detector as a function of jet threshold, 
and corrected for detector effects. The electron and muon measurements were performed in slightly different 
phase space, due to the different rj range and electron-jet isolation requirements, as well as for the different 
QED treatment between electron and muon definitions. The observed signal yields were corrected to recover 
the yield at particle level as described in Section |Vl] 

The corrected ratio i?jct of the production cross sections in the leptonic (electron or muon) decays of the 
gauge bosons W and Z in association with exactly one jet is shown in Figure [3] as a function of the jet px 
threshold for the electron (left) and muon (right) channels. As the jet px threshold increases, the ratio i?jot is 
expected to decrease as the effective scale of the interaction becomes large compared to the difference in boson 
masses. This dependence is observed in the data. The values for the lowest jet px threshold of 30 GeV are: 



i?j^t(e) = 8.73 ± 0.30 (stat) ± 0.40 (syst) 
%t(M) = 8.49 ± 0.23 (stat) ± 0.33 (syst) 

Statistical uncertainties were evaluated by repeating the measurement with Monte Carlo pseudo-experiments 
assuming Poisson distributed data with a mean at the observed yield. Both electron and muon channel results 
are individually compatible with the theoretical predictions. 
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FIG. 3: Results for Rj^t in the electron channel (left) and in the muon channel (right) for their respective fiducial 
regions. The results are compared to NLO predictions from MCFM (corrected to particle level using Pythia). Data 
are shown as black points at the lower bin edge corresponding to the jet px threshold with black error bars indicating 
the statistical uncertainties. The yellow band shows all systematic uncertainties added in quadrature and the green 
band shows statistical and systematic uncertainties added in quadrature. The theory uncertainty (dashed line) shown 
on the MCFM prediction includes uncertainties from PDF and renormalization and factorization scales. Note that these 
threshold data and their associated uncertainties are correlated between bins. 
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Electron and muon channel results were found to be compatible and therefore combined to reduce the statisti- 
cal and uncorrelated systematic uncertainties on the result. Each channel was extrapolated to a common phase 
space, defined as \rip\ < 2.5 before any QED radiation (Born level) with Pythia. The electron channel was 
further corrected for the effect of the electron-jet isolation requirements on the acceptance. This extrapolation 
to a common fiducial region decreases the value of the ratio for both channels primarily due to the more central 
distribution of leptons from the Z. The results were combined using a Bayesian approach jlO) in the combination 
of systematic uncertainties accounting for correlations between them. The systematic uncertainties from E^^^^, 
jet energy scale and resolution and electroweak background sources were considered fully correlated between 
the electron and muon channels. The combined result is shown in Figure |4] (left). The value of iJjot for the 
lowest jet pt threshold of 30 GeVwas found to be 8.29 ± 0.18 (stat) ± 0.28 (syst). This combined measurement 
was also extrapolated to the full phase space, as shown in Figure |4] (right). For a jet px threshold of 30 GeVthis 
ratio was found to be 10.13 ± 0.22 (stat) ± 0.45 (syst). 
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FIG. 4: Left: Combined electron and muon results for Rj^t in a common fiducial region. The results are compared to 
predictions from MCFM (corrected to particle level). Data are shown with black error bars indicating the statistical 
uncertainties. The yellow band shows all systematic uncertainties added in quadrature and the green band shows 
statistical and systematic uncertainties added in quadrature. The theory uncertainty (dashed line) includes contributions 
from PDF and renormalisation and factorization scales. Right: Combined electron and muon results for i?jot extrapolated 
to the total phase space. Note that these threshold data and their associated uncertainties are correlated between bins. 



IX. SUMMARY 



We present a measurement of the ratio of the production cross sections of the gauge bosons W and Z in 
association with exactly one jet, as a function of jet threshold. This ratio was measured in fiducial phase 
space of the detector separately for muons and electrons. These results were also extrapolated to a common 
phase space and combined, as well as extrapolated to be the full phase space of the boson decay products. These 
results were provided as a function of jet pt threshold from 30 to 200 GeV, exploring the transition region of 
electroweak scale breaking in perturbative jet production. The accepted theoretical model was found to be 
consistent with all results. This measurement has the advantage of having reduced theoretical and experimental 
systematic uncertainties due to correlations between the W and Z processes. This measurement builds the 
foundations of a high precision test of the Standard Model, and provides model-independent sensitivity to new 
physics coupling to leptons and jets. Comparisons with LO and NLO perturbative QCD predictions were made 
and found to be in agreement with data over the jet px threshold range covered by this measurement. 
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